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We consider the weak turbulence of whistler waves in the in low- /3 inner magnetosphere of the Earth. Whistler 
waves with frequencies, originating in the ionosphere, propagate radially outward and can trigger nonlinear 
induced scattering by thermal electrons provided the wave energy density is large enough. Nonlinear scattering 
can substantially change the direction of the wave vector of whistler waves and hence the direction of energy 
flux with only a small change in the frequency. A portion of whistler waves return to the ionosphere with a 
smaller perpendicular wave vector resulting in diminished linear damping and enhanced ability to pitch-angle 
scatter trapped electrons. In addition, a portion of the scattered wave packets can be reflected near the 
ionosphere back into the magnetosphere. Through multiple nonlinear scatterings and ionospheric reflections 
a long-lived wave cavity containing turbulent whistler waves can be formed with the appropriate properties 
to efficiently pitch-angle scatter trapped electrons. The primary consequence on the Earths radiation belts is 
to reduce the lifetime of the trapped electron population. 
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I. INTRODUCTION 

The dynamics of whistler wave packets in the inner 
magnetosphere have been investigated using linear ray 
tracing for decades^—. These studies are important in 
understanding the lifetime of trapped energetic electrons 
in the Earth's radiation belts because pitch-angle scatter- 
ing by electromagnetic whistler waves into the loss-cone 
can be the dominant loss mechanism-. The rate of pitch- 
angle scattering is a sensitive function of the spatial, spec- 
tral, and wave-normal distributions of the wave energy 
density, and therefore mechanisms that control the dis- 
tribution of wave energy are important. Recently the re- 
sults of Hasegawa and Chen— and Tripathi, Grebogi, and 
Liu— have been extended to demonstrate that, through 
nonlinear (NL) induced scattering by thermal electrons 
in low-/3 plasmas, quasi-electrostatic lower-hybrid waves 
can scatter into electromagnetic whistler wave s 9 ' 10 . In 
this paper we report on initial efforts to incorporate this 
NL mechanism into the calculation of the propagation 
of whistler waves in the plasmasphere. We find that 
when the turbulent wave energy density exceeds a thresh- 
old, NL scattering can scatter waves that have exhausted 
their usefulness for pitch-angle scattering back into use- 
ful waves capable of further pitch angle scattering. We 
suggest a weak turbulence framework for incorporating 
these effects into ray tracing calculations. This frame- 
work is a natural extension of linear techniques such as 
ray tracing and quasi-linear diffusion that have been the 
workhorse of radiation belt physics for decades. 

Magnetospherically reflecting whistlers (MRW) are 
whistler waves, possibly originating from lightning 
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strikes, whose ray paths in the magnetosphere roughly 
follow magnetic field lines, reflecting many times between 
the magnetic poles. Waves originating in the ionosphere 
propagate to higher and higher altitude and the perpen- 
dicular component of the wave vector increases^. Based 
on linear theory there are two important physical effects 
that limit the degree of influence that MRW are pre- 
dicted to have on the lifetime of trapped energetic elec- 
trons. First, linear damping of whistler waves depends 
strongly on the perpendicular wavelength (increasing as 
kjj. Whistler waves have a significant parallel electric 
field when the perpendicular wavelength becomes compa- 
rable with the electron skin depth. This parallel electric 
field is the driving force of linear Landau damping (res- 
onant interactions with electrons) as well as collisional 
damping. Also, electron drift E x B energy increases as 
the square of k± . Thus as whistler waves propagate from 
their source region and the perpendicular component of 
the wave vector becomes larger the linear damping of 
the waves becomes stronger. Second, the rate of pitch- 
angle scattering of trapped electrons depends strongly 
on the perpendicular wavelength. As the perpendicu- 
lar wavelengths become shorter than the electron skin 
depth the whistler wave becomes more electrostatic in na- 
ture (more lower-hybrid like) and thus the magnetic com- 
ponent (which pitch-angle scatters through the Lorentz 
force, (e/c)v x B) is reduced. In this article we will 
demonstrate that when the energy density of whistler 
waves is large enough, NL induced scattering can return 
a portion of short-wavelength quasi-electrostatic lower- 
hybrid waves back into long- wavelength electromagnetic 
whistler waves which are less damped and more effective 
at pitch-angle scattering. 

The scattered waves (with smaller k±) have only a 
slightly lower frequency than the original whistler waves, 
but the direction of k and their trajectory through the 
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magnetosphere can be dramatically altered*^. A por- 
tion of the scattered waves can be returned back to the 
ionosphere, and some of these waves will retrace the orig- 
inal ray paths back to the magnetosphere. These scat- 
tered waves are then able to efficiently pitch-angle scatter 
rclativistic electrons once again. Through multiple NL 
scatterings in the magnetosphere and reflections in the 
ionosphere a cavity containing turbulent whistler waves 
may be formed that is more long-lived than would be 
predicted without NL scattering, because the effective 
kj_ and thus linear damping rate is reduced. 

In section [TT] we review the dynamics of whistler wave 
packets without NL scattering from the ionosphere into 
the magnetosphere and describe our assumptions in char- 
acterizing these trajectories. In section ITTll we describe 
how we include the effects of NL induced scattering into 
the dynamics of whistler wave packets. In section IIVI 
we describe how a long-lived wave energy cavity may be 
formed due to induced NL scattering and describe some 
of the basic properties of this cavity. In section [V] we dis- 
cuss the role that this newly postulated wave cavity may 
play in radiation belt physics and outline future efforts to 
improve our understanding of the properties of the wave 
cavity. 



II. RAY TRACING WHISTLER WAVES 

A. Model Assumptions 

To model the dynamics of whistler wave packets in the 
magnetosphere/ionosphere we use the following disper- 
sion relation applicable to cold plasmas for waves with 
frequencies between the proton and electron cyclotron 
frequencies^, 
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where, u? LH = *£ iona w 2 s /(l+w 2 e /0 2 ) is the lower-hybrid 
frequency, ui ps is the plasma frequency, and VL S the cy- 
clotron frequency of species s. In Fig. [TJ we show con- 
tours of constant frequency in (k\\,k±,) — (k\\,k±)c/u pe 
space_ calculated from Eq. Q] In the lower right cor- 
ner (k±_ > 1 and ku < k± (me/rrii) 1 ^ 2 ) are the quasi- 
electrostatic lower- hybrid waves with uj 2 ~ w 2 ff , in the 
upper left corner (k± < 1 and k» > {m e /mi) 1 /' 2 ) are the 
electromagnetic whistlers with u> 2 ~ k 2 k 2f l 2 , and in the 
lower left corner (k± < 1 and k\\ < (me/rrii) 1 ^ 2 ) are the 
magnetosonic waves with uj 2 ~ k 2 uj\ H ~ k 2 V\, where 
V\ = B 2 / (Aimimi) is the Alfven velocity. 

We use magnetic coordinates {x, cx, /3) where B = 
Vx = Va x V/3 with a pure dipole field such that, 

Mcos(0) Afsin 2 ((9) „ , 
X = a = (3 = (2) 

where M = —8 x 10 25 G/cm 3 is the earth's dipole mo- 
ment. In these coordinates the components of the group 
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FIG. 1. Contours of constant uj normalized to ujlh calculated 
from Eq. [T] 



velocity may be calculated, 
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where g xx , g aa , and g^P are components of the metric 
tensor. 

To model the density in the magnetosphere we use a 
simple analytical two ion species model that was designed 
to capture the main features of a standard ionospheric 
density profile relevant to whistler wave propagations. 
The densities are all proportional to the magnetic field, 
n s = fs(r)B, multiplied by a function that depends only 
on the radius. The multiplying functions are defined by, 

fe = feo(r - ro) + /ei(l + tanh(^-^)) exp(--^-) 



1/2(1 + tanh( 
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fo = l~f H 



(4) 



with the values of parameters determined by a fit to a 
night-time ionosphere profile. The densities and lower- 
hybrid frequencies are plotted as a function of altitude in 
the upper ionosphere in Fig. [2] 

In addition to following the dynamics of wave-packets 
through (x, k) space, we also follow the amount of energy 
lost due to linear damping. Since the region of interest 
(lower magnetosphere) is a cold plasma (T e ~ 0.1 — 0.5 
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FIG. 2. Densities (top) and lower-hybrid frequency (bottom) 
vs altitude in model ionosphere. 



eV) Landau damping by the thermal part is completely 
negligible compared to collisional damping. There some- 
times exists a population of suprathermal electrons that 
can be in resonance with whistler waves^, however this 
damping is smaller than the collisional damping we con- 
sider (see appendix |A~|) . We include the effects of linear 
collisional (electron-ion) damping by using the approxi- 
mate damping rate, 
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wave vector. The obliqueness of the wave allows for cross- 
field propagation so that wave packets initiated in the 
ionosphere propagate to higher and higher altitude (and 
L-shell). Eventually the wave-packet reaches an L-shell 
greater than the L-shell at which the frequency of the 
wave matches the lower-hybrid frequency and the radial 
group velocity reverses. By this time (a few seconds) the 
wave packet has a large perpendicular wave vector and 
follows the field line more closely. Its radial group veloc- 
ity is relatively small and the wave packet slowly makes 
its way back to the lower-hybrid resonant surface. Fig. 
[3] shows the dynamics of a single wave packet, character- 
istic of packets that reach the magnetosphere, launched 
from an ionospheric source at an altitude of 1100 km at 
30° N and midnight local time with a frequency of 6 kHz 
(where the local lower hybrid frequency is about 12 kHz). 
The model ionosphere/magnetosphere is independent of 
the azimuthal coordinate and thus the azimuthal compo- 
nent of the wave vector is a constant of motion. Therefore 
the wave-packet in figure [3] which was launched with no 
azimuthal component of the wave vector is confined to 
the meridian it began in. In Figure U we show the tra- 
jectory of a wave packet launched from the same point in 
space and the same frequency, but with a component of 
the wave vector out of the noon-midnight plane. In this 
case the packet initially has a significant group velocity 
in the azimuthal direction, but the packet quickly settles 
down into a plane of constant local time and does not 
proceed to encircle the earth. These packets last almost 
10 seconds before the initial energy is reduced by 99% by 
collisional damping. 
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where v e i = 3 x 10~ 5 n e T e 3 ^ 2 (sec~ 1 ), which requires a 
temperature for which we use a simple two temperature 
model profile given by, 

T e = 0.1 + 1/5(1 + tanh( alt ~^ 100 )) eV (6) 

where alt is the altitude in km. 

For this initial work we confine ourselves to consid- 
ering sources in the ionosphere, however, we note that 
geomagnetic storms and chorus emissions are examples 
of potential sources from the outer magnetosphere. As 
for sources in the ionosphere that could generate whistler 
wave power there are many candidates, e.g. lightning, 
VLF transmitters, and velocity ring distributions gen- 
erated naturally or man made such as in the Buaro 
experiment^ or an envisioned future experiment^. 



III. THEORY 

In the framework of weak turbulence theory the dy- 
namics of whistler wave packets in the magnetosphere 
may be described by the wave-kinetic equation which 
gives the evolution of the whistler wave packet "num- 
ber density" N = W/oj where W is the energy density 
and uj is the frequency, 
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~dt 



dr 

dt 



dk 



-N)+V k -^—Nj = lNL N + 2 lL N + Q 

(7) 

where the left side of ([7]) expresses the conservation of 
phase volume in (r, fc) space of wave packets (Liouville 
theorem). The evolution of the wave packets through 
phase space (x, k) are given by the usual ray tracing 
equations, 



B. Linear Dynamics of Whistler Waves 



dr 

dt 



duj dk 
9fc and dt 



dr 



(8) 



Generally speaking, whistler waves with small wave 
normal angles propagate from their source region roughly 
along magnetic field lines, reflecting between the mag- 
netic poles and increasing the value of the perpendicular 



On the right hand side of Eq. (J7J Q is a source of wave 
energy e.g. lightning discharges, chorus waves propagat- 
ing into the plasmasphere, VLF transmitters, etc., 7^ is 
the linear damping/growth rates which along with the 
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FIG. 3. Wave packet trajectory in x-z plane for 6 kHz wave 
launched from 1100 km altiitude and 30° N (top) and evolution 
of k and energy (bottom). 

usual quasi-linear equations for the distribution of parti- 
cles would give the complete self-consistent linear dynam- 
ics of waves and particles. The new physics is contained 
in Inl, the NL induced scattering rate^ilii, 

1 fcf \E kl \ 2 |fei x fc 2 |jj 

Wfe2 1 + fcf ^ iwnm e k^k^n 

1 + fc? Le , ,i I 2 

1 l e fci-fe 2 + £ fei-fe 2 I 

where £ e k '*_ k2 are the plasma electron and ion suscepti- 
bilities, and Eki is related to N k i (detailed discussion is 
contained in the appendix), jnl is to be interpreted as 
the rate of change of the energy contained in waves with 
wave vector (kn,k±) due to all other waves with differ- 
ent frequencies and wave vectors with (fein, fc_u). This 
formula was derived for isothermal, T e ~ T iy low-/3 mag- 
netospheric plasmas from the drift kinetic equation for 
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FIG. 4. Wave packet trajectory for 6 kHz wave launched from 
1100 km altitude and 30°N with component of k in azimuthal 
direction. Top is x-z plane. Bottom is x-y plane. 



electrons with a fluid response for ions. It is the elec- 
tromagnetic generalization of the lower-hybrid induced 
scattering rate found by Hasegawa and Chen— and used 
in tokamak heating scenarios. The scattering rate can 
be physically understood as the landau damping of ther- 
mal electrons by the electric field structure created by 
the beating of two waves. It is important to point out 
that the nature of this electric field structure is inher- 
ently 3D as can be seen by the presence of the cross 
product of two perpendicular wave vectors and the par- 
allel components of the wave vectors within the plasma 
dispersion function. 2D simulations with the simulation 
plane perpendicular to the magnetic field, or 2D simula- 
tions with the magnetic field along one of the simulation 
axes will miss this fast scattering rate^ and will generally 
see only a slower scalar nonlinearity. Three-dimensional 
simulations, or 2D simulations where the magnetic field 
is oblique to the simulation plane are required. For the 
purposes of this paper, the importance of this three- 
dimensionality is that NL scattering allows wave packets 



with components of the wave vector in a meridional plane 
to scatter out of the plane as will be shown in the next 
section. Equation ([9]) was developed in the framework of 
weak turbulence theory, which considers the NL coupling 
of linear waves with stochastic uncorrelated phases, so 
that all waves are solutions to the same linear dispersion 
relation, and hence the NL induced scattering essentially 
redistributes energy in fc-spac o 15 ' 16 . 

From Eq. (|9]) we see that the NL scattering rate is 
largest when the change in frequency Auj^i — Ute|&||i — 
fc|||, and since the electron temperature in the plasmas- 
phere is of the order of an electron volt (which gives an 
electron thermal speed of about 4xl0 7 cm/s) then the 
change in frequency due to a scattering is much smaller 
than the change in wave vector. This allows one to view 
the lines of constant frequency in Fig[T]as all of the possi- 
ble waves that a given wave could scatter into. For exam- 
ple, if we have a wave that is lower-hybrid like, marked 
with point (a) in the figure, it may scatter into a whistler 
wave, marked with point (b) in the figure. This has im- 
portant consequences for the influence of whistler waves 
on the radiation belts, namely that waves that have al- 
ready become lower-hybrid like due to the evolution of in- 
dividual wave packets in a dipole field as predicted by ray 
tracing, have an opportunity to become whistler waves 
again. 

To gain a detailed understanding of the effects of the 
NL scattering on the propagation of wave power in the 
radiation belts a self-consistent solution, either numeri- 
cal or analytical, to the wave-kinetic equation is needed. 
Since this is beyond the scope of this work, we adopt 
a test wave-packet approach to investigate the effects of 
NL scattering and thus solutions to Eq. [71 



IV. EFFECT OF NL SCATTERING ON WAVE PACKET 
TRAJECTORIES 

In this section we describe how a cavity in the magne- 
tosphere may be formed when the wave energy density 
exceeds a threshold, we describe the properties of the 
cavity, and estimate the influence the formation of the 
cavity will have on the trajectory and fate of whistler 
wave packets and on the lifetime of trapped energetic 
electrons. 



A. Access to magnetospheric cavity from ionospheric 
source 

We choose a source region in the ionosphere or in the 
transition region between the ionosphere and the mag- 
netosphere (see Table fl} and suppose that the whistler 
energy flux is nearly monochromatic and isotropic in k 
space. Launching whistler wave packets with the same 
frequency, from the same location in space, with random 
directions of the wave vector we can build up a back- 
ground. By choosing a random direction for the wave 
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FIG. 5. Trajectories of wave packets that enter magneto- 
spheric cavity launched from 1100 km altitude, 30° N, with a 
frequency of 6 kHz. 

vector many packets will immediately be lost by propaga- 
tion towards the Earth's surface, and another percentage 
of packets will be lost by closely following the magnetic 
field lines and reaching the Earth's surface in the oppo- 
site hemisphere. In addition, some of the packets are 
trapped by the steep electron density gradient (closest 
to the peak in density near 400 km), and because of the 
high density are quickly collisionally damped. We call 
these packets ducted. The percentages of packets for the 
different cases are listed in table HI In Fig. [5] we show 
the trajectories of wave packets that enter the magne- 
tospheric cavity launched from 1100 km altitude, 30° N 
with a frequency of 6 kHz. This figure gives an idea of 
the three dimensional geometry of the cavity which will 
be further discussed in the section IIVCI In Fig. [5] we 
show a sample of wave packet trajectories launched from 
500 km altitude with a frequency of 6 kHz that are lost 
to the Earth. We find that the percentage of packets that 
enter the magnetospheric cavity increases for the higher 
altitude and latitude initial conditions, mainly because 
more packets that travel to the southern hemisphere are 
refracted back into the magnetospheric cavity. 

B. Threshold for NL scattering 

We estimate the threshold in wave energy density nec- 
essary for NL induced scattering to play an important 
role in the evolution of whistler waves. To do this we 
suppose that the source lasts for a time much longer than 
the bounce time of wave packets between the magnetic 
poles and follow wave packets until they reach the surface 
of the earth or their energy density reduces by 99% while 
oscillating near the lower hybrid surface in the magne- 
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FIG. 6. Trajectories of wave packets that are lost to the Earth 
launched from 500 km altitude, 30° N, with a frequency of 6 
kHz. Red paths are lost to the Northern hemisphere. Black 
paths are lost to the Southern hemisphere. 



tosphere. From this background of waves we can form a 
grid in space and calculate ensemble averages of useful 
quantities such as k±. We dchnc the following average, 

y wx 

(X) ( Xi ) = "g^ w (10) 

where W is the energy of the particular wave-packet. In 
Figure [7] we plot such quantities from an ensemble of 
waves launched from 500 km altitude with a frequency of 
4 kHz (63% of the LH frequency) . Since we assume that 
the source lasts for 30 seconds we record all equatorial 
plane crossings at different times to form this average. 

Next, we estimate the NL scattering rate (Eq. [9]) by as- 
suming a broad band of turbulence, Suj ~ cj, and replace 
the summation with integrations. Using the condition 
that the integrand is dominated by the condition £ ~ 1, 
and the change in frequency is comparable to the ion- 
cyclotron frequency (see Appendix ) the NL scattering 
rate may be approximated as, 




Then, from the form of Eq. (O we expect that NL scat- 
tering will dominate linear damping when the energy den- 
sity is large enough such that ^nl > 1l- Performing 
the same averaging procedure on the collisional damping 
rate, Eq. ([5]), we can estimate the NL threshold as, 

(Bi\ 1 v ei /mx3/2 / l + (fc^) 2 \ 5uBl 
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FIG. 7. Ensemble average of the perpendicular wave vector 
for waves launched from 500 km altitude with a frequency of 
4 kHz. 

in the limit (kn ) < which is justified by the calcula- 
tions leading to Fig. [7J Then taking the peak (k±) from 
the ray tracing calculations, and estimating Suj/uj ~ 1/2 
we can then estimate the NL threshold energy. We find 
that the energy density to reach threshold corresponds 
to magnetic field amplitudes of about 100 pT. It should 
be pointed out that this estimated threshold is of the 
same order of magnitude as the observed amplitude of 
plasmaspheric his s 17 ' 18 . Thus it can play an important 
role in such phenomenon as radiation belt slot formation 
and it can effect theories of the source of plasmaspheric 
hiss. As we will show in the next section NL scattering 
can significantly effect the wave propagation. 



C. Properties of Wave Energy Cavity 

In the previous sections we computed the accessiblity 
of whistler waves into the magnetosphere and the thresh- 
old in local wave energy density necessary for NL scatter- 
ing to be important. Now we use ray tracing calculations 
to estimate the volume of the cavity so that we may es- 
timate the energy necessary to be input to the cavity to 
achieve sufficient NL scattering. We launch 1000 wave 
packets into the magnetosphere and compute the maxi- 
mum extent in the azimuthal direction (in units of hours, 
i.e. 12 hours is 180 degrees), and the maximum extent 
in latitude. In table [TT] we summarize the results of our 
calculations. In general we find that the azimuthal ex- 
tent of the cavity decreases with increasing frequency and 
altitude, and is fairly insensitive to latitude. The merid- 
ional extent of the cavity increases with the latitude and 
altitude of the whistler launch point and frequency. To 
compute the maximum radial extent we find the the max- 
imum radial value achieved by at least 10% of the packets. 
The reason for this distinction is that for the higher lat- 
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FIG. 8. Wave packet trajectory (top) for 4 kHz wave launched 
from 750 km altiitude and 30°N (blue) with wave vector con- 
fined to meridional plane and scattered after 1.4 seconds (red). 
Perpendicular wave vector and energy of wave as a function 
of time (bottom), blue indicates unscattered wave and red is 
scattered wave. 



itude cases there are a minority of packets (< 10% that 
can achieve large radial excursions on the first or second 
magnetospheric reflection). We find that in general the 
maximum radial extent decreases with frequency and in- 
creases with latitude. Also included in the table is an 
estimate of the volume of the cavity where the minimum 
radius of the cavity is always taken to be 1.4 Re. From 
these volume estimates we can then use the wave energy 
density necessary to overcome the NL threshold and ar- 
rive at the necessary total energy to achieve threshold. 
This is indicated in the final column of table HU 

Next we describe an example of a single scattering. For 
this purpose we assume that enough wave energy density 
is present in the system so that a wave packet may scat- 
ter. In Fig. [5] we show the trajectory (top in blue) of 



a wave packet with a frequency of 4 kHz launched from 
750 km altitude at 30° N latitude. This wave packet has 
the wave vector in the noon-midnight meridian so that 
the wave packet is confined to the same meridian. After 
1.4 seconds the wave packet has increased its value of k± 
to about unity and experiences a large angle scattering, 
such that the new wave packet has a component of the 
wave vector out of the meridian plane. The new scattered 
packet also has k± ~ 1, but this packet is on a different 
trajectory so that k± decreases in time for a while be- 
fore returning back to the magnetospheric lower hybrid 
surface and increasing k± again. This reduction in k± 
allows for the energy density in the wave to be almost 
a factor of 2 larger than the energy density found in the 
packet without scattering (and a subsequent reduction in 
the final value of k±). 

Finally, we discuss the effects of many scatterings. In 
figure [8] we show what happens to a particular wave 
packet that is scattered which says in the magnetospheric 
cavity. However not all packets stay in the cavity. Some 
are lost to the Earth, and others experience large linear 
damping. To begin to quantify this energy loss we de- 
fine a whislter wave albedo that we calculate by launch- 
ing whistler wave packets with random directions for the 
wave vectors from inside the magnetospheric cavity and 
integrating the ray tracing equations for a specified time 
and then comparing the total energy of the packets at 
this time to the total initial energy in the packets. Thus, 
the albedo is defined as, 

Albedo = ^ t = t ' nd , (13) 

where the numerator does not include packets that are 
lost to the earth. A good choice for the ending time 
in the albedo calculation would be the time a particular 
packet spends before it's next "collision", i.e. 1/jnl- 
For simplicity, we chose the energy density such that the 
NL collision time was 1.5 seconds for all packets. In Fig. 
IHlwe show a sample of ray tracing calculations for which 
the wave packets remain in the cavity, with their initial 
conditions marked by circles and their termination points 
by x's. Performing these kinds of calculations for all of 
the cases we considered we summarize the results in table 
HID In general we find that the albedo decreases with 
increasing frequency, is fairly insensitive to the altitude of 
release, and is larger for the higher latitude release case. 
In addition to the albedo the table shows the percentage 
of packets that are lost to the Earth. In all cases this is 
the dominant source of energy loss. 

Another feature of many NL scatterings is that the 
size of the cavity can increase, this feature can be seen 
in the example figure [9j where the size of the new cavity 
can be seen by the ray trajectories in blue overlaid on 
the original size of the cavity estimated as thick red line 
in the figure. One mechanism for increasing the size is 
that packets which have become confined to a particular 
meridian may be scattered such that they propagate az- 
imuthally. An individual case can be seen in Fig. [5J In 
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a » 1, the pitch-angle diffusion coefficient becomes, 



A if 



W 



(B 2 /(8tt) 



cos(6»)cot 3 (6i) (15) 



which can be very small. What this formula shows, is 
that without NL scattering of whistlers the time with 
which a wave packet can effectively scatter relativistic 
electrons is very short. However, with NL scattering a 
population of packets can be scattered back so that the 
increase in k± from ray tracing is reversed, and the rate 
of pitch-angle scattering can remain significant. 



FIG. 9. Scattered wave packet trajectories for 6 kHz wave 
launched from 1100 km altiitude and 30° N. Solid thick lines 
(red) indicate approximate boundary of unscattered cavity. 



table Mil we estimate the dimensions of the cavity after a 
first scattering. 



D. Effect of Cavity Formation on Energetic Trapped 
Electrons 



There are three main components to enhancing the 
pitch-angle scattering of energetic trapped electrons by 
NL scattering of whistler waves. The first is that the 
energy density of the whistler waves from a given source 
can be enhanced above the value predicted by ray tracing 
without NL scattering. This was demonstrated in figure 
IH1 a corollary of this phenomenon is that the waves can 
live longer in the cavity with NL scattering than without, 
and this was discussed in relation to the albedo of the 
Earth. The final ingredient is that the average wave- 
normal angle can be decreased, which enhances the pitch- 
angle scattering in the following way. 

The pitch-angle diffusion coefficient is approximately 
given by, 



Dn 



W 



J! (a) 



(14) 



where a — k^v±/(fl e /"fn) and is the relativistic fac- 
tor. The cyclotron resonance condition {uj — k\\V\\ = 
f2 e /7i?) for relativistic electrons and whistler waves gives 
fc||ccos(0) ~ Qe/jR where 6 is the electron pitch-angle. 
With this condition the argument of the Bessel function 
becomes a ~ (k±/ku) tan(#) which quickly becomes large 
due to linear ray tracing of wave packets. In this limit, 



V. CONCLUSIONS 

In conclusion, we have shown that NL induced scatter- 
ing can alter the propagation of whistler wave energy in 
the magnetosphere from ionospheric sources, such that 
a long-lived cavity may form that can more effectively 
pitch-angle scatter trapped relativistic electrons and af- 
fect the lifetime of the trapped population. We have 
introduced a framework for incorporating this NL effect 
into future more detailed studies, and have used ray trac- 
ing calculations to explore the properties of the NL cav- 
ity. Next we highlight a couple of observations from ra- 
diation belt physics for which NL scattering may play an 
important role. 

Plasmaspheric hiss, which plays an important role 
in pitch-angle scattering, refers to observed broadband, 
electromagnetic, incoherent, and turbulent fluctuations 
on the whistler branch mostly confined to the high- 
density, cold plasmasphere. Typical amplitudes of the 
magnetic fluctuations range from 10 to 100 pTiL±£, i.e. 
in the range where the NL threshold could be met. The 
observation of the wave normal angles of hiss is the most 
contentious and difficult to make. Wave normal angles 
have been observed to be (1) rather oblique, within 45 
degrees of perpendicular to the magnetic field 19 , (2) par- 
allel propagating 4 (in both directions), and (3) a mixture 
of both types (1) and (2 ) 20 ' 21 . Small wave-normal angles 
of plasmaspheric hiss have been difficult to explain theo- 
retically, because ray tracing studies (as shown in Section 
IIIBI) show that k± becomes large very quickly. The most 
established mechanism involves refraction of whistlers off 
of the density gradient at the plasmapause, however there 
are observations where the plasmapause may be very 
distant. The presence of these small-wave normal an- 
gle whistlers may be due to NL scattering provided the 
energy density is sufficiently high. The source of plasma- 
spheric hiss has been attributed to lightnin g 22 i 23 , chorus 
emissions outside of the plasmaspher o 24 ' 25 , and amplifi- 
cation of background or embryonic waves by anisotropic 
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Altitude 


Latitude 


Frequency 


LH Frequency 


% into Cavity 


% into North 


% into South 


% Ducted 


(km) 


(° N) 


(kHz) 


(kHz) 


(%) 


(%) 


(%) 


(%) 


500 


30 


4 


6.3 


8.5 


63 


26 


2 


500 


30 


6 


6.3 


8.0 


57 


35 


0.4 


1100 


30 


4 


12 


41 


52 


7.6 





1100 


30 


6 


12 


38 


52 


10 





1100 


30 


8 


12 


38 


52 


10 





750 


30 


7 


8.4 


14 


52 


35 





850 


30 


7 


9.7 


18 


51 


31 





950 


30 


7 


10.9 


25 


52 


23 





750 


45 


7 


9.9 


47 


53 








500 


45 


4 


7.5 


35 


63 





3 


750 


45 


8 


9.9 


46 


54 








600 


45 


6 


8.2 


48 


52 









TABLE I. Access to magnetospheric cavity computed from 1,000 ray tracing calculations for each case. 



Altitude 


Latitude 


Frequency 


Azimuthal Extent 


Maximum Radius 


Latitudinal Extent 


Volume 


Energy 


(km) 


n 


(kHz) 


(Hrs) 


{Re) 


(°) 


(R%) 


(MJ) 


500 


30 


4 


2.5 


2.3 


61 


2.1 


0.5 


500 


30 


6 


2.0 


2.2 


65 


1.4 


0.4 


1100 


30 


4 


2.3 


2.4 


64 


2.4 


0.6 


1100 


30 


6 


1.9 


2.2 


68 


1.6 


0.4 


1100 


30 


8 


1.7 


2.2 


71 


1.3 


0.3 


750 


30 


7 


1.9 


2.2 


68 


1.3 


0.3 


850 


30 


7 


1.9 


2.2 


69 


1.5 


0.4 


950 


30 


7 


1.8 


2.2 


69 


1.4 


0.4 


750 


45 


7 


2.2 


3.1 


89 


7.4 


1.9 


500 


45 


4 


2.6 


3.3 


88 


10.8 


2.7 


750 


45 


8 


2.2 


3.0 


89 


6.9 


1.7 


600 


45 


6 


2.4 


3.0 


89 


7.0 


1.8 



TABLE II. Dimension of magnetospheric cavity computed from 1,000 ray tracing calculations for each case 



electron distribution s 17 ! 26 , and the debate about which 
source is more likely has hinged in some part to the az- 
imuthal accessibility and the wave-normal angle of waves 
all computed from linear ray tracing. As we have demon- 
strated both the wave-normal angle and the azimuthal 
extent from a given source can be affected by NL scat- 
tering. 

As evidence of the importance of the wave-normal an- 
gle of plasmaspheric whistler waves on the lifetime of en- 
ergetic trapped electrons Meredith et ai— found that 
the lifetime of 2-6 MeV electrons measured by SAM- 
PEX in the inner slot region, that was filled due to 
the Halloween storms of 2003, was barely influenced by 
un-ducted whistler waves (with predicted large wave- 
normal angles), however by including a small popula- 
tion of whistler waves with small wave normal angles (as- 
sumed to be ducted whistlers) the observed lifetimes of 
these electrons could be explained. 



In the future, we plan to develop a more self-consistent 
solution to the wave kinetic equation to investigate the 
evolution and equilibrium of spectral wave-power in the 
radiation belts. To compare these to observational data 
more detailed models of the sources and equilibrium 
structure of plasma density, and magnetic field would 
be useful. Finally, we mention one future research topic 
which could have profound consequences on the ampli- 
tude of turbulent whistler waves and that is the inclusion 
of the growth of whistlers due to anisotropic electron dis- 
tributions (loss-cone distributions) 28 by the mechanism 
put forward in application to magnetospheric physics by 
Kennel and Petschek— . The criterion for growth by 
this mechanism has been difficult to identify because, as 
was demonstrated in Section III B[ whistler waves quickly 
develop short perpendicular wavelengths and experience 
large linear damping, and thus waves would have only 
a short time to gain energy from the unstable electron 
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Altitude 


Latitude 


Frequency 


Lost 


Albedo 


Azimuthal 


Maximum 


Latitudinal 


Volume 












Extent 


Radius 


Extent 




(km) 


n 


(kHz) 


% 


% 


(Hrs) 


(Re) 


n 


(R%) 


500 


30 


4 


22 


65 


4.6 


2.8 


93.5 


10.6 


500 


30 


6 


47 


38 


3.5 


2.5 


88 


5.4 


1100 


30 


4 


19 


68 


4.3 


2.9 


99 


12 


1100 


30 


6 


39 


45 


3.4 


2.6 


95 


6.4 


1100 


30 


8 


57 


28 


2.6 


2.4 


88 


3.5 


750 


30 


7 


51 


34 


3.1 


2.5 


87 


4.4 


850 


30 


7 


50 


34 


3.1 


2.5 


89 


4.5 


950 


30 


7 


50 


34 


3.4 


2.5 


89 


5 


750 


45 


7 


32 


49 


4 


3.0 


115 


13.4 


500 


45 


4 


17 


70 


5.4 


3.0 


122 


21 


750 


45 


8 


37 


43 


3.2 


2.9 


110 


9.4 


600 


45 


6 


29 


53 


4 


3 


112 


14 



TABLE III. Properties of magnetic cavity after the first scattering. 



population. However, with NL scattering whistlers can 
spend a longer time with properties desirable for growth 
and thus may better be able to tap into the unstable elec- 
tron population. Detailed analysis of this possibility will 
be done in the future. 



Appendix A: Landau and Transit-Time Damping due to 
Suprathermal Electrons 

If Landau damping due to the suprathermal electrons 
(100 eV to about 1.5 KeV) is sufficiently large, then the 
wave energy required to reach the threshold for NL scat- 
tering can potentially be very large. In this appendix 
we calculate the collisionless (Landau and transit-time) 
damping caused by the population of suprathermal elec- 
trons and show that it is smaller than the collisional 
damping we have considered. Consequently, for the pa- 
rameters of interest the suprathermal electrons will not 
play any significant role in the formation or lifetime of 
the cavity within the plasmasphere. 

Suprathermal electrons are observed as a sporadic 
phenomenon-^. To evaluate the damping rates we use 
properties of the observed statistical average population 
of suprathermal electrons. Most studies focus on data 
from higher-L shells than what are examined in our work 
(most of the ray-tracing calculations in our study are con- 
fined to L<2.5). For example, Thorne and Horned model 
the suprathermal population with a series of Maxwellians 
between 4.2<L<5.7. Bell et al^ model the suprather- 
mal electron population between 2.3<L<5.7 with power- 
laws and find that the density of suprathermal electrons 
is about 200 times less (at 150 eV) than at the higher 
L-shells. Recently Li et al£^ reported on global distri- 
butions of suprathermal electrons and confirmed that for 
energies less than a few KeV the results of Bell et al. are 
more appropriate for inside the plasmapause and that the 



density increases with increasing L-shell. Since our stud- 
ies are at L-shells below these observational works, it is 
likely that our calculations of collisionless damping over- 
estimate the damping rates. Finally, it should be pointed 
out that if there is a population of suprathermal electrons 
contributing to Landau damping, then consequently the 
lifetime of these electrons in the radiation belts would be 
short as these particles would be energized into the loss 
cone and precipitated into the ionosphere. Because of the 
sporadic nature of suprathermal flux there will be periods 
with no resonant electrons to cause collisionless damping. 
However, to be conservative in the following we assume 
that the population is always present and contributing 
to collisionless damping 

We calculate the damping rate by first calculating the 
wave power dissipated by the energetic but tenous pop- 
ulation of suprathermal electrons, 

P = E-J=^-E*- Xa -E (Al) 

where x a ^ s the anti-hermitian part of the susceptibility 
due to the suprathermal electrons. By considering only 
Landau damping (n = 0), taking 7 < oj and k±p e << 1 
then the limit of the susceptibility tensor given by Stix 32 
can be used to calculate the power dissipated as 

Hnuj J \ fc|| J dv\\ 

(A2) 

where the first term in the bracket is transit-time damp- 
ing (from /iV-B mirror force), the last term is Landau 
damping due to parallel electric field, and the middle 
term is the combination of these terms. Next we stipu- 
late that the suprathermal electrons (with a ratio of at 
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most n s th/n 10~ 5 ) do not significantly alter the fields as 
found using the cold plasma approximation and the lim- 
its used to derive the dispersion relation of Eq. (). The 
fields may then be related as, 



kuk 



lr-L 



-ifc 2 



1 



k\ 



UJ 



kj_c 



E, E, E,, E„ — - — B z 

. (A3) 

where uj = u;/|fi e | and k = kc/uj pe . The power dissipated 
by the suprathermal electrons is then entirely expressed 
in terms of the parallel component of the magnetic field, 



P = -U)- 



\B, 



Sn 



4W£ 



2irv±dv± 



v ± v \\ 
V Ae 



k i co 



dvnS I <" 



k 2 k 4 



UJ 

h 



dfo 

dvn 



kl) 



UJ 



(A4) 



where Vke = c\ft e \/uj pe ~ 5 x 10 9 cm/s. For Maxwellian 
suprathermal distributions these integrals can easily be 
done, yielding, 



P 



I 2 o vl 



v: 



Air 1 / 2 k 



1 - 



<h (1 



ii Ktn (i + fci)' 



■ exp 



Ae 



J sth 



v Ae 



2k 2 



< th (1 



4fc 4 



(A5) 



where typical values of the ratio VAe/vsth ~ 3 — 10. By 
examination of Eq. (|A5I) the transit time damping is typ- 
ically a little smaller than Landau damping, and interest- 
ingly acts in a way as to reduce the damping rate. Bell 
et al^ have published an analytical distribution function 
given by, 



ft 



Bell 



a 

7a 



c 

,,6 



(A6) 



where a, b, and c are constants. The integrals in Eq. 
(|A7|) can also be calculated with Eq. (|A6[) . assuming the 
functions go to zero more quickly than 1 /w 4 outside of the 
range where data has been collected. We have computed 
these integrals and present the numerical values of the 
damping rates resulting from these calculations below in 
table El 

For the Landau and transit time damping to be signif- 
icant the phase velocity of the waves must coincide with 
the velocity of the suprathermal electrons. This limits 
wave- vectors to specific ranges to consider. For example, 
whistler waves (with (me/mi) 1 / 2 < ku < 1 and k± < 1) 
have phase velocities that match electrons with energies 
between 100 and 1500 ev. In this case, the dispersion re- 
lation simplifies to Co 2 = k 2 k 2 . Then to get the damping 
rate 7 = — P/W we need the energy density of the waves, 



W 



1 

167T 



B * -B + E*^-e h E 

ouj 



(A7) 



which for whistler waves is dominated by the mag- 
netic field energy, becoming simply W ~ (I-BjJ 2 + 
|£ z |7(16tt) ~ (1 + fc(j/A:i)|£ 2 |7(87r). Then the damp- 
ing rate of whistler waves with [ra e jmi) x / 2 < ku ~ fcj_ < 
ldue to Maxwellian suprathermal electrons becomes, 



risth Vl 



Ae ;„5 



exp 



J sth 



V 2 

,,2 



V 



si h 



(A8) 



where the terms proportional to v 2 th / (2k 2 V] Le ) < 1 have 
been dropped. On the other hand, quasi-electrostatic 
lower-hybrid waves (k± < 1 and feii < fci(m e /m 1 ) 1 '' 2 ) 
have phase velocities that match the velocity of electrons 
with 10's of KeV, for which there aren't many suprather- 
mal electrons. Thus, toward the magnetospheric edge of 
the cavity, where such waves are prevalent, the collisional 
damping will always dominate. 

Finally, we quantitatively compare the collisionlcss 
damping rates to the collisional damping rate for different 
values of (ku,k±) in Table HVl In this table we calculate 
the collisionless damping rate in two ways. JL—Bell is 
calculated from Eq. (|A4[) . where we have computed the 
velocity space integrals using Bell's published analytical 
distribution function given by Eq. (|A6[) . and dividing by 
Eq. (IA7|) where we have used the cold plasma expres- 
sions for th and the fields in Eq. (|A3|) . jl- Maxwellian is 
calculated from Eq. (| A5|) and Eq. (|A7I) . For whistlers 
this rate becomes Eq. (|A8[> . For a suprathermal density 
we took n s th = 0.013 and v s th corresponding to 150 eV. 
For each of these calculations we choose typical L=2 pa- 
rameters of B = 0.04 G and n — 5 x 10 3 /cm 3 . The first 
case (fc||,fcj_) = (1/7, 1/7) corresponds to whistlers inside 
the cavity, where we see that the collisionless damping 
is at least one order of magnitude smaller than the colli- 
sional damping we have considered. The second case cor- 
responds to a magnetosonic wave, and the third and the 
fourth cases correspond to a lower-hybrid like wave char- 
acteristic of the magnetospheric edge of the cavity. In 
these last three cases the damping rate due to suprather- 
mal electrons is at least an order of magnitude smaller 
than the collisional damping we have considered. 



Appendix B: Estimation of NL Scattering Rate 

In this appendix we simplify the NL induced scattering 
rate of the main text. We start with a general form of 
the NL induced scattering rate as published in Ganguli 
et al.— , 
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|fei x k 2 
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k\ 



(Bl) 
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fclj 


kx 


JColl (rad/s) 


E Res (KeV) 


-JL-Beii (rad/s) 


iL-(Maxwellian) (rad/s) 


1/7 


1/7 


3xl0" 2 


0.3 


5xl0" 5 


8xl0~ 4 


0.02 


0.1 


4xl0~ 3 


0.2 


lxlO" 5 


3xl0~ 5 


1/7 


1/2 


lxlO" 1 


1.4 


2xl0~ 3 


2xl0" 4 


1/7 


1.5 


0.3 


1.9 


2xl0~ 2 


5xl0~ 5 



TABLE IV. Damping rates assuming B = 0.04 G and n = 5 x 10 3 /cm 3 . For 

"/L-(Maxweiiian) we take n st h = 0.013/cc and V Bt h 

corresponding to 150 eV. 



where, assuming a fluid model for the ions and drift ki- 
netic model for the electrons, 



-fci-fca 



-fci-fca 



E 



ft 2 _ ( Aw )2 



(B2) 



(Afc) 2 ti/2 



(i + C^K)) 



and is a real function while e| _ fe has an imag- 

inary part from the dispersion function. Here Aw = 
— Wfc 2 , Afe = |fei — k 2 \ and Ak z — fc lz — fc 2z are 
the frequency and wave vectors of the beat wave, and 
the subscripts 1 and 2 denote the mother and daughter 
waves respectively. For NL wave scattering by magne- 
tized electrons the condition 



c 



I k 2z ~ fa J 



< 1 



(B3) 



must be satisfied, i.e. there must be enough particles that 
resonate with the beat wave structure and C-^(C) ~ 1. To 
maximize the NL rate the inequality \e k _ k \ > |e| _ k | 
is necessary. This inequality is satisfied under two cases: 



(1) Q s < |Aw| < AkC s -> {Akfpl > 1 



(B4) 



is the subsonic condition (meaning the beat wave has a 
frequency below the ion-sound wave frequency), and 



(2)|Aw-ni| = {Akfp 2 s 



(Akfp 2 < 1 (B5) 



is the ion cyclotron condition (meaning the beat wave has 
a frequency below the electrostatic ion cyclotron wave 
frequency), where C 2 = 2T e /m,i and p s = C s /fli. When 
(|B3[) and one of the conditions (|B4[) or (|B5[) is satisfied 
the scattering becomes, 



fi 2 



1NL 



k 2 



\E kl \ 2 |fei x fe 2 ||| 



U!k2 1 + ^2 t 



^2 47mT e 



P i- 2 



l + fc 2 



(B6) 



Considering a broadband spectra of width 5w 3> Aw we 
simplify the sum over all k by considering only the waves 
contributing to the summation given the conditions of 
Eq. (|B3]) and (|B4jt or (|B5]> , Since the NL scattering rate 
is maximum for large angle scatterings we estimate that 
\ki x k 2 \ 2 / (k^k^) ~ 1- For the subsonic case (case 1 



above) the NL scattering rate may be estimated by the 
fraction of interacting waves (Aw/<5w ~ AkC s / Suj) under 
the condition of Eq. jB3) and pi]). 



1NL 



where. 
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W fc2 1 + ^2 n T e 1 + tfaY 
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Afc/3, 



1/2 



(B7) 



nT e 



\E kl \ 
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ft 2 ^ 8nnT e 87tnT ( 

e ki 



and ( fci ) is the average value of fa . Note that Nuj = W\ 
and Eq. (|B8|) gives the relation between the plasmon 
number density and the amplitude of the fluctuating 
magnetic field. 

For the ion cyclotron case (case 2 above) the beat wave 
frequency is the ion-cyclotron frequency, but the fei's 
that contribute to the summation are restricted by the 
condition (|B5I) . Afc2p2 < i ; so that we can write the NL 
scattering rate as 



f2 2 k\ 



1NL TTT2 m Ak 

w fc2 1 + fef n T e i + /fay dw \m e 

(B9) 

Now, we evaluate these conditions near the lower-hybrid 
resonant surface in the magnetosphere to determine 
whether the rate for scattering of lower-hybrid waves to 
magnetosonic/whistler waves is governed by Eq. (|BT|) 
or (|B9p . Here the plasma is almost all hydrogen and 
f) e ~ 6 x 10~ 5 so that (Ak) 2 p 2 sI1 = (Ak) 2 (3 e m H /m e ~ 
O.l(Afc) 2 . In this case, to meet the condition of Eq. (|B4[) 
would require (Ak) 2 > 10, which from ray-tracing (as 
can be seen in Figure [7]) cannot be satisfied . However it 
is possible to meet the condition (|B5|) since (Ak) 2 < 10, 
and the condition (Eq. 



Aw 



fi 



H 



\k 2z ~k lz \v te Ak z v te 



< 1 



(B10) 



which evaluates to (Afc z ) 2 > 10m e /mH, is not too re- 
strictive. Therefore Eq. (|B9|) is the proper estimate 
of the NL scattering rate for the formation of the mag- 
netospheric cavity, which is somewhat smaller than the 
rate as would be given by Eq. (|B7I) . The rate given by 
Eq. (|B7I) is applicable to the scattering of lower hybrid 
waves to magnetosonic/whistler waves in the ionosphere 
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where oxygen plasma dominates and the main parameter 
Pe m o/ m e is of the order of unity. 

To arrive at Eq. (fTTj) in the text we consider the large- 
angle NL scattering from lower-hybrid waves to whistler 
waves so that Ak ~ (k) in Eq. (|B9|) . From ray-tracing 
waves with u) ~ ujlh we find that kj_\ > fciii and we 
estimate the scattered wave vector fc_i_2 ~ (k±i). With 
these assumptions and Eq. (|B8I) we arrive at 



m I 



3/2 



5a; 



l + (fcxi) 2 ) 



2 R 2 
^0 



(Bll) 



where the subscript 1 has been dropped in Eq. (|11|) . 
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